clear;clc;
close all;
FimageQ=load('Fimage_Qacoustic.mat');

FimageA=load( 'Fimage_acoustic.mat');

FimageA_Qdata=load( 'Fimage_acoustic_attenuated.mat');
[zz,xx]=size(FimageQ.Fimage);
dx=8;
dz=8;
x_tick=[0:50:xx];x_label=dx*x_tick;
z_tick=[0:50:zz];z_label=dz*z_tick;
Font=14;

plotimage(FimageA.Fimage);title('acoustic');
set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.6]);
axis([0,xx,0,zz]);set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x (m)','FontSize',Font);
ylabel('z (m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);

plotimage(FimageA_Qdata.Fimage);title('acoustic with attenuated data');
set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.6]);
axis([0,xx,0,zz]);set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x (m)','FontSize',Font);
ylabel('z (m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);

plotimage(-FimageQ.Fimage);title('Q acoustic');
set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.6]);
axis([0,xx,0,zz]);set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('x (m)','FontSize',Font);
ylabel('z (m)','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'XTick', x_tick);
set(gca,'XTickLabel',x_label,'FontSize',Font);
set(gca, 'YTick', z_tick,'FontSize',Font);
set(gca,'YTickLabel',z_label);


%%%%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%%%
%%%            plot trace             %%%%%%%%%%
%%%%^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^%%%
xi=150;
figure;set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.4]);
plot(FimageA.Fimage(:,xi),'r','LineWidth',1);
hold on;
plot(-FimageQ.Fimage(:,xi)*0.5,'--b','LineWidth',1);
plot(FimageA_Qdata.Fimage(:,xi),'k','LineWidth',1);
xlim([0,zz]);%set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('z (m)','FontSize',Font);
ylabel('Amplitude','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'xTick', z_tick,'FontSize',Font);
set(gca,'xTickLabel',z_label);

xi=200;
figure;set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.4]);
plot(FimageA.Fimage(:,xi),'r','LineWidth',1);
hold on;
plot(-FimageQ.Fimage(:,xi)*0.5,'--b','LineWidth',1);
plot(FimageA_Qdata.Fimage(:,xi),'k','LineWidth',1);
xlim([0,zz]);%set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('z (m)','FontSize',Font);
ylabel('Amplitude','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'xTick', z_tick,'FontSize',Font);
set(gca,'xTickLabel',z_label);

xi=250;
figure;set(gcf,'units','normalized','position',[0.1 0.1 0.6 0.4]);
plot(FimageA.Fimage(:,xi),'r','LineWidth',1);
hold on;
plot(-FimageQ.Fimage(:,xi)*0.5,'--b','LineWidth',1);
plot(FimageA_Qdata.Fimage(:,xi),'k','LineWidth',1);
xlim([0,zz]);%set(gca,'xaxislocation','top');
set(gca,'tickdir','out');
xlabel('z (m)','FontSize',Font);
ylabel('Amplitude','FontSize',Font);
set(gca,'FontSize',Font);
set(gca, 'xTick', z_tick,'FontSize',Font);
set(gca,'xTickLabel',z_label);
